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We observe critical phenomena in spherically symmetric gravitational collapse of 
Einstein Cluster. We show analytically that the collapse evolution ends either in 
formation of a black hole or in dispersal depending on the values of initial parameters 
which characterize initial density and angular momentum of the collapsing cloud. 
Near the threshold of black hole formation, we obtain scaling relation for the mass 
of the black hole and find the critical exponent value to be 3/2. We numerically 
confirm that there exist wide ranges of initial parameter values around the critical 
configuration for which the model remains shell-crossing free. 



I. INTRODUCTION 

Critical behavior in gravitational collapse was first found by Choptuik in numerical sim- 
ulations of spherically symmetric massless scalar field [l| . For marginally supercritical data, 
it was found that the black hole mass scales as M ~ — p*)'^ , where p is a parameter of 
family of the initial data which characterizes strength of the configuration, and p* is the 
critical value. The critical exponent 7 has a value 0.37 for scalar field, and is universal in 
the sense that it is the same for all one parameter families of initial data. 

Subsequently, the matter models were generalized. Critical behavior was found in the 
collapse of radiation fluid [2I, perfect fluid with p = kp 3], and in the collapse simulations 
of Yang Mills field with critical exponent 0.19 [J] and also in axion-dilaton model from low 
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energy effective string theory Although critical behavior has been found in numerical 
simulations of many matter models, very few analytical examples are available so far. Koike, 
Hara and Adachi gave explanation of the scaling of black hole mass and the universality by 



carrying out renormalization group studies These studies made it clear that the critical 
exponent for a given model can be obtained by finding the growth rate of the unique unstable 
mode of the self similar critical solution which is the fixed point of the renormalization 
group transformation. Subsequently, linear stability analysis became a standard technique 
to demonstrate critical behavior and calculate critical exponent in gravitational collapse 
models. On the other hand, it is not necessarily clear where and how the linear order 
eigenvalue analysis gets applicable for initial data highly nonlinearly deviated from the 
critical solution, prior to numerical simulations. 

Exact value of critical exponent has been calculated analytically in few models, which 
include a matter model with thin shell coupled with null fluid [7|, massless scalar field in four 
dimensions js] and in (2+1) Ads Spacetime [ol. We present here a model in which threshold 
of black hole formation can be seen and corresponding exponent can be calculated exactly 



by a different method 10|. 



Here we consider a spherically symmetric collapsing model in which non-colliding particles 
move in such a way that the radial pressure vanishes and non-zero tangential pressure is 



present in the collapsing cloud. Such a static system was first introduced by Einstein 



which was later generalized to non-static case Il2l. In recent years, the properties of this 

nn 

model have been studied extensively p^, I14I |. 

We observe that the model shows criticality near the threshold of black hole formation. 
We work in the comoving coordinates and investigate evolutions of individual shells. It is 
analytically shown that the collapsing cloud either forms a black hole or disperses depending 
on the values of initial parameters which characterize the initial density and tangential 
pressure of the collapsing cloud. We derive scaling relation for the mass of the black hole 
near its threshold. 

This model has a limitation that the radial stress is vanishing and the Misner-Sharp 
mass for the collapsing cloud is time independent. However, this simplifies the Einstein 
equations and make the model tractable. The model considered here is different from the 
models considered for the criticality so far in the following aspects. Critical exponent is 
obtained without invoking self similarity. We have found the critical self similar solution for 
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this system and it will be discussed elsewhere |15[|. The no-shell crossing conditions make 
the initial data restrictive, and only certain ranges of initial data is allowed in the comoving 
coordinate system. The model shows critical behavior for two parameters as long as we 
assume regularity at the initial epoch and the smoothness of the initial parameters, and 
thus, in the space of initial data sets, the requirement of black hole formation within this 
model already restricts the space of possible sets to be of codimension one, in contrast with 
the standard critical behavior where the shape of the initial data have very large variety. 
Despite all these limitations, the model remains interesting as one can see transparently how 
an infinitesimal mass black hole formation can take place when we fine tune the initial data. 

The outline of the paper is as follows. In Section II, we discuss collapse equations and 
regularity conditions and in Section III, the Einstein cluster model is discussed. In Section 
IV, critical behavior is investigated and it is demonstrated how a given sets of initial value 
parameters decide the final outcome. In Section V scaling law for the black hole mass is 
derived. Shell crossing is discussed in Section VI and conclusions are outlined in Section 
VII. 

II. EINSTEIN EQUATIONS, REGULARITY AND ENERGY CONDITIONS 

We use the polar coordinates (t, r, 6', 0) to write the spherically symmetric metric as 

ds' = -e^'^i'^^'Ut' + e2^(*'^)rfr2 + R^t, r)d^l\ (1) 

where dVl' is the line element on two-sphere. We also take the energy-momentum tensor to 
be diagonal for the collapsing Type /matter field (that is, the frame is a comoving coordinate 
system) which is given by 

Tl = -p- T;=p.; T^ = T^ = po. (2) 
This is a fairly general class of matter fields, which includes various known physical forms of 



matter 16| . The quantities p, Pr and pe are density, radial pressure and tangential pressure 
respectively. We take the matter field to satisfy the weak energy condition, that is, the 
energy density as measured by any local observer be non-negative, and so for any timelike 
vector we have 

TikV'V" > 0. (3) 
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This amounts to 

p > 0; p + p, > 0; p + pe > 0. (4) 

The dynamical evolution of the system is determined by the Einstein equations, and for the 
metric ([T]) in the units SttG = c = 1, these are given as 

F' F 

R' 

jy'{p + Pr) ='2{Pe-Pr)^ -P'r, (6) 

- 2R' + R'^ + ^§ = 0, (7) 
G-H = l-^, (8) 

where (') and (') represent partial derivative with respect to t and r respectively and 

G{t,r) = e-^^{R')\ H{t,r) = e'^^R^ (9) 
The quantity F{t,r)/2 is the Misner-Sharp mass for the collapsing cloud, which gives total 



mass within a shell of comoving radius r at time t [17|. In order to preserve the regularity 
at the initial epoch, F{ti, 0) = 0, that is, the mass function should vanish at the center of 
the cloud. It can be seen from the equation ([5]) that density of the matter blows up when 
R = or R' = 0. Here the case R' = corresponds to the shell- crossing singularities. 
Now let us write physical radius as 

R{t,r) = rv{t,r). (10) 

Using the scaling independence of the coordinate r and initial collapse condition, we write 

v{U, r) = 1; v{ts{r), r) = 0; v{U, r) < 0, (11) 

where ti and ts stand for the initial and the singular epochs respectively. The condition 
v{ti,r) < signifies that we are dealing with initially collapsing shells. We scale the radial 
coordinate r in such a way that at the initial epoch R = r, and at the singularity, R = 0. 
The advantage of the introduction of this new variable v is that the regular center at r = 
(where we also have i? = 0) is now distinguished from the genuine singularity at i? = 0. We 
now have w = 1 at the initial epoch, and f = at the singular epoch R = 0, but at all other 
epochs in-between v has a non-zero finite value for all values of r. 



5 



III. NON-STATIC EINSTEIN CLUSTER MODEL 



The spherically symmetric collapse models, where the radial pressure is taken to be 

pressure could be non-zero have been studied in quite some 



vanishing and the tangentia 
detail over past many years 13|, UM, llSl]. 

The Einstein cluster is an example of such a cloud where tangential stresses are present. 
This is a spherically symmetric cluster of rotating particles where the motion of the particles 
is sustained by an angular momentum which has an average effect of creating a non-zero 
tangential stress within the cloud. Neighboring shell particles are counter-rotating such that 
spherical symmetry is preserved. 

We consider such a non-static cluster of gravitating particles in four dimensions. For the 



non-static Einstein cluster models the equation of state is given by 

1 / L2 



12| 



2\R^ + L\'P^ 

where L(r) is a function of the radial coordinate r only and is known as specific angular mo- 
mentum. Vanishing radial pressure implies that the Misner-Sharp mass is time independent. 
Regularity of the initial density at the center requires 

F{r)=r^M{r), (13) 

where M.{r) is a smooth function. It is clear that as ^ 0, p — oo. Thus the density blows 
up at the singularity i? = which will be a curvature singularity as expected. Let us now 
define a suitably different iable function A{r^v) in the following manner 

v\r,v) = A{r,v)^^B!. (14) 
Then from equation ([6]), we have the equation of state given as 

pe = ^A^vRp- (15) 
Now using the equation ([H]), we can integrate ([7]) to get 

G = b{r)e^'^. (16) 

Here b{r) is another arbitrary function of the comoving coordinate r. Following a comparison 
with dust collapse models we can write 

b{r) = l + r\{r), (17) 
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where bQ{r) is the energy distribution function for the collapsing shells. Finally, using equa- 
tions (fll, ([161) and (flTl) in (El) we have 



A comparison of equation (|T2l) with (|T5l) gives 

^''^ " R{R^ + L^y ^^^^ 
We can integrate above equation to get 

Considering initial density, pressure and energy profiles to be smooth would ensure L(r) also 
to be smooth. 

We write the initial profiles in the form 

M{r) = mo + m2r^ H , (21) 

L^(r)=iy + kr^ + ---, (22) 
bo{r) = boo + bo2r^ + ■■■ ■ (23) 

The initial density profile, initial specific angular momentum profile and velocity profile of 
the cloud are chosen fully in equations ( I2T1) . ( !22l) and ( !23l) respectively. Regularity at the 
initial epoch requires leading order term of to go as r^. Now we need to evolve this initial 
data according to equation (ITSl) . 



IV. CRITICAL PHENOMENA 

Initially, at the onset of gravitational collapse all the shells have the scale factor v{ti,r) 
value as unity and v(ti, r) < 0, which implies an initially collapsing cloud. Bounce of a shell 
is indicated by the change in the sign of v. Evolution of a particular shell may be deduced 
from equation (fT8|) . Rewriting (ITSl) in terms of the scale factor we have 

V(r, v) = (^v^ + vv\ (24) 

We call V{r, v) as the effective potential for a shell. It can be expressed in terms of initial 
profile functions of the system as follows 
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V{r, v) = -[hoV^ + Mv^-^v+ '-^ ] . (25) 



The allowed regions of motion correspond to V{r, v) < 0, as vv"^ is non-negative. The first 
factor in equation (12^ is always positive, because it is the g^^ term of the metric tensor and 
the quantity in parenthesis in the same equation is also always positive, hence location of 
the turning points (where i) changes its sign) are not decided by these two terms. The main 
features of the evolution of a shell basically derive from the cubic polynomial in equation 

The dynamics of the shells may be studied by finding the concerned turning points. If 
we start from an initially collapsing state {v < 0), we will have rebounce for a shell if we 
get V = 0, before the shell has become singular. This can happen when V{r,v) = (see 
for example Fig.l). Hence, to study the various evolutions for a particular shell we must 
analyze the roots of the equation V{r,v) = 0, keeping the value r to be fixed. The method 
employed here is similar to another class of tangential pressure collapse models where the 



io|. 



metric function qqq is taken to be a function of the physical radius R 

Out of three roots of the cubic polynomial only positive real roots correspond to 
physical cases. We take 6o to be positive. The region between the unique positive 
roots is forbidden since in those regions v'^ < 0. For a particular shell to bounce, it 
must therefore lie, during initial epoch (f = 1), in a region to the right of the second pos- 
itive root. We can see that three different types of evolution of a particular shell are possible: 



i) If the effective potential for a shell has two distinct positive roots in the range [0, 1], 
then the shell bounces off. 

ii) If V{r,v) < in the whole range [0,1], then the shell will reach the singularity at 
v = 0. 

iii) When the potential has a double root in [0, 1], then it indicates that the shell is in 
critical collapse condition. 



Let us assume that we are working in initial data space for which there are no shell- 
crossings. We will discuss the issue of shell- crossings in the matter cloud in §6. If a particular 
shell with comoving coordinate bounces then all the shells with coordinate r > must 
also bounce. This implies that to investigate the situation when the entire cloud is just 
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about to disperse off, it is sufficient to study the dispersal of the shells near the center. 
Therefore, to find the threshold of black hole formation and to get the scahng relation for 
the mass of the black hole, we need to analyze the model only near the center. We have 
neglected here higher order terms in r in the expansion, since we want to consider only the 
evolution of the shells near r = 0. 

To analyze the dynamics in this model in detail, we study one particular configuration 
of the initial data, in which mo, m2 are kept fixed, 600 is taken positive, Z4 is set to zero and 
Iq is allowed to vary. Non-zero value of Z4 will make the central shell and hence the cloud 
to always bounce and black hole will not be formed. Close to the center, we neglect higher 
order terms in the expansion, in which case effective potential can be written as 

V ~ — 6o(^) ~ {iTT'O + Tn2r^) + l^r^ v — (mo + m2r^)l6r^, (26) 
where bo{r) — boo + bo2r^- We use an arbitrary unit since there is no length scale in the 
original system of equations. For a cubic equation ax^ + bx^ + ac + d = 0, a double root 
occurs when discriminant A for the polynomial vanishes, where 

A = ld>abcd - Ab^d + 6^c^ - Aac^ - Tla^d^. (27) 

For above cubic equation this takes the form 

A = r^Ai(r,p?^), (28) 

where are the initial parameters. Now we study collapse evolutions for three different 
categories of the initial data and p\. 

A) Supercritical Evolution 

Fig.l shows the effective potential for the central shell (r = 0) and for three outer shells 
(r = 0.10,0.113,0.15) for a small value of /g. Potential for the central shell has always a 
double root (r = 0, Ai(0,p|) 7^ 0). If we go slightly away from the center, potential for a 
shell goes negative for the entire range of v which allows the shell to reach the singularity 

If we increase r further, the potential maxima starts going up, and there is a value of r 
at which the effective potential again has a double root. We call this radius as the critical 



9 



radius of the collapsing cloud for the chosen set of the initial numbers. At this value the 
discriminant of the cubic vanishes (Ai(rc,p|) = 0). 
Apparent horizon is given by 

F = R, (29) 

which near the center in (f,r) plane can be written as Vah{r) = r'^mo, where at f = Vah a 
shell becomes trapped. It can be seen from the figure that if the potential remains negative 
in the range [0,1], the value of v decreases and finally goes below Vah which indicates the 
trapping of the shell. Now as the shells below Vc go inside the apparent horizon, a black 
hole eventually forms. These shells finally reach the singularity and contribute for the mass 
of the black hole formed, but the shells with comoving coordinate more than critical radius 
bounce off. As the black hole forms for these values of Iq, this represents supercritical region. 




FIG. 1: Figure represents effective potential for different comoving radius in supercritical configu- 
ration. Here 1^ = 3.5, bo = 4.0, mo = 1.15 and m2 = —20. 



B) Critical Evolution 



As we increase /g further, we reach a value at which the potential for the central shell 
has double root (r = 0,Ai(0,p^) = 0) and all other shells have two distinct real positive 
roots. We call this value as the critical value l^c of the initial parameter /g for the other 
fixed parameters. This situation is shown in Fig.2. The outer shells have positive potential 
and forbidden region, therefore, all those shells will bounce back. This configuration gives 
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the critical solution of the system. The critical point is the boundary point between the 
dispersal and the situation when the collapse to form black hole just begins. 
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FIG. 2: Figure represents effective potential for different comoving radius in critical configuration. 
Here Iq = 5.3, 60 = 4.0, mo = 1.15 and 1712 = —20. 

C) Subcritical Evolution 

If we increase Iq further to l^c potential goes even more positive. There is a complete 
bounce of the collapsing shells. Hence this initial configuration generates subcritical evo- 
lution. However it should be noted that only the central shell reaches the singularity 
(Ai(0,p^) 7^ 0) and all other shells bounce. This massless central singularity is a timelike 
naked singularity which can be seen as follows. For the central shell e'''^^'^^vv'^ = b^QV + uiq, 
and as mo and 600 are positive, central shell reach singularity at f = while all outer 
shells bounce off at value of v which is larger than Vah for those shells and are not trapped. 
Thus, central singularity forms but trapped surfaces do not form in the cloud, making the 
singularity globally visible 

V. MASS SCALING LAW 



r=0 

r=0.10 

r=0.15 





To derive the scaling law for black hole masses, it is the key to estimate the critical shell 
radius, which separates the completely collapsing cloud from the surrounding dispersive 
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FIG. 3: Figure represents effective potential for different comoving radius in subcritical configura- 
tion. It can be seen that outer shells bounce. Here /g = 10, 6o = 4.0, mo = 1.15 and m2 = —20. 

cloud. The expression for the critical radius in terms of the initial parameters can be obtained 
in the supercritical region from the condition that at critical radius, effective potential 
has two equal roots. The double root condition for polynomial in fl26l) together with the 
consideration of shells only near the center gives 



ArriQ 



ml Iq 



-16m^m2 - 18 6oo"^o^6 + 2 mom2k + 4 fcoo^e 



(30) 



At = he = 4mQ, the critical radius vanishes. It can be seen in Fig.4 that as we increase le 
in the supercritical region, the critical radius monotonically decreases from a positive value 
and becomes zero at /gc- If we increase Iq further all the non-central shells in the cloud 
bounce. Therefore, we can write 



2 

7^ 



k \Ir — I 



6c 



(31) 



where A; is a constant. In the tangential pressure model, the Misner-Sharp mass depends 
only on r. Mass which collapsed to form singularity from the regular initial profile is given 

by 



(32) 
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FIG. 4: Critical radius in the supercritical region. The curve is the locus of the point on which Ai 
vanishes. 

As we are analyzing the system near the center and also near the threshold of black hole 
formation, critical radius ^ 0, therefore the mass of such black hole can be written as 

Mbh ^ Cm\k - kcf^ , (33) 

where Cm is a constant near the threshold. 

We can fix any of the parameters mo and and vary the other one to obtain the expression 
for critical radius. Following the same procedure as depicted earlier, it can be easily seen 
that for mo, the same scaling relation and critical exponent exists near the threshold of the 
black hole formation. Therefore, in general for the parameter rj, we write 

Mbh ^ Cr,\r] - r]c\^ . (34) 

VI. SHELL CROSSINGS 

Above analysis is based on the condition that for the chosen initial data near the critical 
values, there are no shell crossings during the evolution of the collapsing cloud, which implies 
that if a shell bounces then all the shells with a larger value of comoving radius will also 
bounce. It is important to check that there are indeed no shell crossings taking place for 
a given initial data. We know that if R'{r,t) > there are no shell crossings, therefore, 
we need to obtain R' for a given set of initial conditions. Equations (|T^ and (fT8|) are two 
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coupled equations for v and R{r,t). They can not be solved explicitly analytically in the 
coordinate system used and numerically also solving them is not easy. 

However, there is a way out, Gair has given a coordinate transformation which eliminates 



the function u in the evo 



ution equation, using which one can obtain R' at constant time t 



in the following manner 19|. He has used a time coordinate r which he calls proper time 



experienced by a dust particle 



After changing the coordinates, we get from equation (ITS]) 

We can choose a common origin of time for all shells, i.e. at some time, t = 0, we set r = 
for all r (see section 3 in [19]). This gives us the initial conditions to integrate above equation 
numerically. We note that 



l-Sm^f^l. (37) 



3r(^)/(^),0/ 

Now, to get the desired quantity what remains is to know {dT/dr)t. To obtain it, we first 
differentiate equation fl35l) with respect to r, then change the order of partial derivatives and 
use equation (l37ll . obtaining a differential equation for {dT/dr)t which is as the following 



where ^ is ^(^^^^^i-^ and ip is ^ (^1 + ^ j . The initial condition for the above equation is 
provided by the fact that we have set r = at t = for all the shells, which implies 
{dT/dr)t = at r = 0. We integrate fl36l) first and then substitute the values of {dR/dr)r 
and {dR/dT)r in equation (15^ . after integrating which equation (155]) gives {dR/dr)t. 



We find that for the chosen initial data in the vicinity of critical situation, there are no 
shell crossings. Physical radius R and {dR/dr)t for the subcritical and supercritical case are 
plotted in the following figures. We do not give here a general criteria or condition on the 
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FIG. 5: Physical radius of different shells in supercritical evolution, where initial parameter values 
are mo = 1.15, bo = 1.95, 62 = 50, mi = 5.0 and L = 5.0. It can be seen that shells with comoving 
radius less than the critical radius reach singularity while those with larger than critical radius 
bounce. 
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FIG. 6: {dR/dr)t in supercritical evolution, where the initial parameter values are mo = 1.15, bo = 
1.95, 62 = 50, mi = 5.0 and L = 5.0. It can be seen that the value of R' at constant time t remains 
positive. 



initial parameter functions which will ensure no shell crossings. However, one can choose 
the initial data which gives no shell crossings and one can see that the critical solution lies 
in this allowed range of parameters. 
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FIG. 7: Physical radius of different shells in subcritical evolution, where initial parameter values 
are mo = 1.15, 60 = 1-95, 62 = 50, mi = 5.0 and L = 5.7. It can be seen that all the shells bounce. 
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FIG. 8: {dR/dr)t in subcritical evolution, where the initial parameter values are mp = 1.15, 60 = 
1.95, 62 = 50, mi = 5.0 and L = 5.7. It can be seen that the value of R' at constant time t remains 
positive. 

VII. CONCLUSIONS 

In this paper, we have presented a simple analytic model of gravitational dynamics which 
is physically motivated in the sense that it is a special realization of the Einstein- Vlasov 
(collisionless particles) system. The system satisfies standard requirements of physical rea- 
sonableness, namely, it has an equation of state, and satisfies energy and regularity condi- 

















r.0.1 


















r.0.2 




- 














r=0.3 


/- 


- 
















- 












































































16 



tions. 

We have shown that both black hole and dispersal are possible outcomes in the gravita- 
tional collapse of this model. The mass of the black hole near the threshold of its formation 
shows power law behavior with critical exponent 3/2. Very small mass black holes could be 
formed if the initial data is fine tuned. The model shows " universal" behavior with respect to 
the parameters which fix the initial density and specific angular momentum of the cloud. As 
shown in the plots, for a generic class of the regular initial data involving central singularity 
formation, shell crossing does not occur. The effective potential method for calculating the 
critical exponent applies to this well as to an another class of tangential pressure 

models jlOl]. This suggests that the method applies at least for all mass-conserving systems. 
However, the system considered here has a limitation and bit different from the models con- 
sidered so far for critical behavior in one respect that it has no radial stress which enables 
shells to interact directly with each other. However, the critical behavior seen in general 
relativity appears to be a generic phenomenon and it is expected that only the radial stress 
is not essential for the phenomena to occur. It will be very interesting to explore whether 
the same method holds for studying critical behavior in the model with radial pressure as 
well. 
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